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Abstract 

•i-H . 

^ , A large fraction of top quark events in pp collisions at 1.8 TeV will contain 

additional soft hadronic jets from gluon bremsstrahlung off the quarks and glu- 
ons in the hard processes qq,gg — > tt bbW + W~. These extra jets can cause 
complications when attempting to reconstruct mj from the invariant mass of 
combinations of final-state quarks and leptons. We show how such soft radia- 
tion cannot be unambiguously associated with either initial-state radiation or 
or with final-state radiation off the b quarks. The top quarks can radiate too, 
and in fact the pattern of radiation has a very rich structure, which depends on 
the orientation of the final-state particles with respect to each other and with 
respect to the beam. We calculate the full radiation pattern of soft jets in the 
soft gluon approximation and compare with several approximate forms which 
are characteristic of parton shower Monte Carlos. The implications for top mass 
measurements are discussed. 



1 Introduction 



Very recently the CDF collaboration has reported evidence for a top quark of mass 
174±17 GeV/c 2 in pp collisions at 1.8 TeV [Q. The leading-order production processes 
are qq, gg — > tt — > bbW + W~ , and a statistically significant signal is seen in the channels 
W + W~ — * Zz//z/ and W + W~ lu + jets (7 = e,/i). From the latter sample, the 
top mass is reconstructed from the final-state lepton and jet momenta. A potential 
problem with this procedure arises when the final state contains an additional hadronic 
jet, since it is not known a priori whether this jet arises from initial-state radiation, 
in which case it should be ignored, or whether it originated (say) in bremsstrahlung 
off a b quark, in which case it should be included in the m t reconstruction. We note 
that a substantial fraction of the candidate top events reported by CDF do indeed 
have an additional jet over and above the number associated with the leading-order 
production and decay process. 

In fact there is no meaningful separation of such jets into 'initial-' and 'final-state' 
radiation ||. In any part of phase space, the soft jet cross section contains contri- 
butions from radiation off the incoming qq,gg, off the produced tt at the production 
stage, off the t and b quarks in the t — > bW decay process, and from interferences 
between these. In specific regions, however, we might expect certain types of radi- 
ation to be more important than others. For example, we would expect radiation 
close to the beam axis to be dominated by gluon emission off the incoming partons, 
while radiation close to a b quark should be dominated by collinear gluon emission off 
that quark. But away from these special collinear regions, all particles (including the 
top quarks) emit gluons with approximately equal strength. In what follows, we will 
study the distribution of soft gluon radiation over all of phase space, to see to what 
extent our expectations about initial- and final-state radiation are valid. We will see 
that while the initial-/final-state radiation picture is too naive, a decomposition of the 
radiation pattern into contributions associated with tt production, decay, and their 
interference provides the information necessary for top mass reconstruction. 

Our purpose in this paper is twofold. First, we wish to study the correct distribu- 
tionsQ to determine where the gluons come from and where they go, in a way that is 
directly relevant to top momentum reconstruction and mass measurement. Our second 
aim is to compare these results with simpler models which are characteristic of the way 
gluon radiation is implemented in Monte Carlo event generators, such as the initial- 
/final-state picture mentioned above. The difference between these approximations 
and the exact treatment could then provide an indication of the errors implicit in 
using certain Monte Carlos to generate extra soft jets in tt production. 

1 Our calculations are 'exact' in the sense that we include gluon radiation off all colored particles 
in the production and decay process. Our only approximation is to assume that the gluon is 'soft', 
i.e. Eg/E <C 1, where E ~ m t is a typical energy scale of the subprocess. 
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2 General formalism for soft gluon radiation 



2.1 Gluon radiation in tt production and decay 

We start with the leading-order process 

a(h) + b(k 2 ) - t( qi ) + i(q 2 ) - b( Pl ) + W + + b(p 2 ) + W~ , (1) 

where ab = qq or gg, and the particles' momenta are indicated in parentheses. Naively, 
the invariant mass of each b — W system is equal to m t . Now suppose that the final 
state contains an extra gluon, with momentum In terms of Feynman diagrams, this 
gluon can be emitted off (taking the qq process as an example) either of the incoming 
light quarks, the s-channel gluon, off either of the top quarks before they decay weakly 
(i.e., on the timescale of the strong production process), off either of the top quarks 
on the timescale of their weak decay, or off either of the final-state b quarks (and also 
from the decay products of a hadronically decaying W). All of these amplitudes can 
of course interfere. 

Because of the infra-red divergence associated with gluon emission, an extra jet in 
a ti event will usually be soft. In this case, one can analyze the radiation pattern in 
the soft gluon approximation, as was done in Refs. H |J. One can write 

1 da a s „ _ 

da dE g d cos 9 g d(p g An 2 9 

where da^ is the differential cross section for the lowest-order process (i.e., with no 
gluon radiation), E g is the energy of the soft gluon, and a s is the strong coupling. The 
function T is the sum of 'antenna patterns' of the radiation from the different sources 
listed above. It can be written generically as 

•F = 3~ PROD 4" J~ DEC ~t~ J~] [NT (3) 

where J^prod is the contribution from emission at the tt production stage (including 
initial-state radiation), ,F DE c is the contribution from emission off the t and t and their 
decay products at the weak decay stage, and F mT refers to the interferences between 
these emissions. More specifically, .F INT contains contributions from the interference 
(i) between radiation in the t decay and that in the i decay and (ii) between radiation 
in ti production and radiation in either decay.Q Note that for soft gluons with E g ~ 
r t ~ 1 GeV, JF INT is sensitive to the top decay width, as discussed in Ref. 0. However, 
the observable soft jets that are relevant to the pp collider experiments have energies 
much larger than T t , and so in practice the interference terms are numerically small. 

2 In terms of the results of Ref. (3], ^prod corresponds to \A\ 2 , T OEC to \B\\ 2 + \B 2 \ 2 , and T mT 
to -2Rc[BiB^] + 2Re[A(B 2 - B x )*]. 
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Explicit expressions for T for the qq and gg subprocesses have been presented in 
and are listed again here for completeness: 



/'prod 

J- DEC 
J-INT 



+ c 2 k 1 q 1 + c z k x q 2 + c 4 fc 2 <?i + c 5 k 2 q 2 + c 6 q^q 2 + c n q x q x + c 8 q 2 q 2 , 
C7[qm + VYP\ - 2glpi] + c ? \q 2 q 2 + p 2 p 2 - 2q^p 2 ] , 

Xi {c^hpi - kxqx) + c 4 [k 2 pi - k 2 q x ] + c & [q^pi - q^q 2 ] + 2c 7 [qtf) 1 - gTgi]} 
+ X2 {c 3 [kip 2 - hq 2 ] + c 5 [k 2 p 2 - k 2 q 2 ] + c 6 [q^p 2 - q^q 2 ] + 2c s [qtf) 2 - q^q 2 ]j 

+ Xl2 C 6 [P1P2 - QlP2 - Q2Pl + gi?2] , (4) 

where the antennae pq and 'profile functions' Xi are defined by 

p-q 



pq 



Xi 



Xl2 



1,2) 



p ■ k q ■ k ' 

( qi -ky + m^ t 1 

m 2 T 2 (qi ■ k q 2 ■ k + rn 2 T 2 ) 

[( qi ■ kf + m 2 v 2 } i(q 2 ■ ky + m?r?] 



(5) 



In terms of the gluon energy, we have JF ~ -^g" 2 ? an d so the cross section has the 
infra-red behavior da /dE g ~ Eg 1 , as expected. Additional collinear singularities arise 
from the p ■ k denominators when p 2 



0. 



The antenna coefficients q depend on the color structure of the hard scattering 
and are different for the qq and gg processes (C F = 4/3, iV = 3): 



coefficient 


— > tt 


#5- -> 


Cl 

c 2 
c 3 
c 4 

c 5 
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c 7 
c 8 


1 

AT 
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2 

N 
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N 

2CV - i 
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V 
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-2C F + 2N + 2Y 
C F -X-Y 
C F + X-Y 
C F + X -Y 
C F -X-Y 
2Y 
-Cf 
-Cf 



The quantities X and K depend on the gg — > ti subprocess energy and scattering 
angle. Explicit expressions and a discussion can be found in [|J. 

In general terms, the antennae which contribute to JFprod correspond to color 
strings stretched between the initial-state partons and the final-state t quarks before 
they decay, and those which contribute to JFdec correspond to color strings linking the 
final-state t{t) and b(b) quarks at the weak decay stage. We remind the reader that 
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color antennae exhibit the 'string effect': radiation between the two momenta in an 
antenna is enhanced compared to that outside them; see for example Ref . || . We will 
see several instances below where this effect influences gluon distributions. 

For purposes of top momentum reconstruction and mass measurement, this decom- 
position serves as a guide to what to do with additional soft jets from gluons. Gluons 
from the production stage are emitted before the t and t quarks go on shell and should 
not be included in top momentum reconstruction. Those from the decay stage are to 
be included along with the top decay products. For the interference contributions no 
such clear-cut assignment is possible; however for the examples we consider here, the 
interference contribution is negligible. 

2.2 Gluon radiation in simpler models 

There are several simpler models which we could construct if we were interested in 
simulating gluon radiation in top events without using the full distribution, Eq. (|). 
Consider, for example, a model which includes only initial-state radiation. This would 
be appropriate, for example, for the purely electroweak process qq — ► Z° — ► l + l~, 
where the color string is stretched between the incoming quarks. In the language of 
the antennae introduced above, this corresponds to 

T = 2C F (2N) k^k 2 , (6) 

for qq(gg) — > tt respectively. We will use T8R' to refer to this model. 

A more sophisticated approximation would also include final-state radiation by 
allowing the b quarks to radiate, assuming that they were linked by a second color 
string, as for example for the process qq — > Z° — > bb. In this case we would findQ 

T = 2C F {2N) hk 2 + C F \2<m>2-<M>x ~KPi\ ■ (7) 

We will refer to this model as 'ISR/FSR'. The ISR and ISR/FSR models are easily 
implemented and are characteristic of what appears in parton-shower Monte Carlos. 

A third possible model corresponds to the naive expectation that the top quarks 
do not radiate because of their short lifetime. (In our case, where the relevant gluon 
energies are larger than r t , the top quarks do radiate; see [[J for a discussion.) This 
model corresponds to taking T t — > oo (and hence Xii X12 — >• 1) in Eq. (|5|) and gives 

T = c x k x k 2 + c 2 k 1 p 1 + c 3 hp 2 + c A k 2 pi + c 5 k 2 p 2 + c 6 prp 2 + c 7 p'rp 1 + c 8 p^p 2 . (8) 

This is the radiation pattern we would see if the b and b were produced directly, and 
so we will refer to this model as 'BB'. 

3 This would in fact be the correct result for the process qq(gg) — > H — > tt — * bbW + W~ in the 
limit T t -> 00 j§. 
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In what follows, we will compare the complete distribution given in Eq. (|]) with 
those from the ISR, ISR/FSR, and BB models, corresponding to Eqs. (0), and 



3 Numerical results and discussion 

In the previous section we have seen how the radiation pattern depends on the ori- 
entation of the final-state t and b quarks. This obviously varies from event to event, 
and in our calculations we must therefore integrate over the full phase space for tt 
production and decay and weight by the appropriate parton distributions. This will 
then yield the correct radiation pattern in the laboratory frame. In particular, we 
will see important differences between the models discussed at the end of the previ- 
ous section. To understand how and why these differences arise, it is helpful to first 
work in the parton subprocess center-of-mass, fixing the t and b momenta in 'typical' 
configurations, as was done in Refs. || |J. 

In order to assess the typical final-state configurations, we first generate a sample of 
events corresponding to qq,gg — > tt — > bbW + W~ — > bbllvv with m t = 174 GeV/c 2 in 
pp collisions at 1.8 TeV. Figure 1 shows the resulting distributions in b(b) transverse 
momentum and pseudorapidity, and in the b — b azimuthal angle difference. The 
majority of b quarks are produced centrally (i.e., within \r}b\ < 1) and with p b T ~ 
30 — 100 GeV. The distribution in azimuthal angle difference shows a slight preference 
for back-to-back production in the transverse plane. We note for future reference that 
the centrality of the fe's is due in large part to the fact that they receive a boost 
in the direction of their parent top quarks, which tend to be produced centrally; 
this also accounts for their slight back-to-back preference. We note also that for 
m t = 174 GeV/c 2 , the cross section is dominated by the qq — > tt process, which is an 
order of magnitude larger than the gg — > tt contribution. 



3.1 Features of the gluon distribution in the tt center of mass 

To investigate the radiation pattern of soft gluon emission, we choose a configuration 
with the following properties. For the moment we work in the tt center-of-mass frame, 
with subprocess energy vl = 3m t . The t and t are produced centrally (j] t — r\i — 0) 
and the b quarks have rjb = —t% and A(p b i = 180°, with the direction of the b defining 
= 0°. For purposes of illustration, we will consider the particular case rjf, — 1, which 
corresponds to p\ = p\ = 57 GeV/c. The distribution of gluon radiation in the rj — <p 
plane is then obtained from Eq. (|2|): 

dN 11 da a, E„ T 
= — i i . (9) 

drjg d(f) g cosh 2 r] g dao dE g d cos 9 g d(f) g An 2 cosh 2 r] g 
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At this stage we are only concerned with the general features of the radiation pattern, 
and so we (somewhat arbitrarily) set E g = 10 GeV and a s = 0.1. 

Figure 2 shows the distribution of Eq. (Q), for the choice T]b = 1. The two most 
obvious features are: (i) the strong enhancement of the radiation near the b and b 
directions, and (ii) the approximately constant 'pedestal' distribution away from the b 
quark directions and in particular at large positive and negative r] g . Closer inspection 
reveals another important feature - the radiation has minima close to the original 
t-quark directions, in this case at r) g ~ 0. These represent the 'dead cones' of the t and 
i quarks 0. In order to study the various individual contributions to the distribution, 
and to compare different approximations, it is more useful to consider particular slices 
through the two-dimensional plot. In Fig. 3 we show the same distribution as a 
function of r) g for the particular values <ft g = 0°, 45°, 90°. The curves correspond to 
the decomposition of Eq. (|J): the long-dashed line is the 'production' contribution, 
the short-dashed line is the 'decay' contribution, the dotted line is the interference 
contribution, and the solid line is the total. For 4> g {= 4>b = 4>t) — 0°, we clearly see the 
dead cone of the b quark (in the decay contribution) and the much broader dead cone 
of the t quark (in the production contribution). The interference contribution is small. 
Away from the 6-quark direction, the production contribution dominates. At <p g = 45° 
the effects of the final-state quarks on the radiation pattern are much reduced, and 
at <p g = 90° - which is exactly half-way between the b and b quarks - the production 
contribution dominates and the distribution is approximately uniform. 

Figure 4 shows the total contribution of Fig. 3 (solid line) compared to the dis- 
tributions corresponding to the ISR (dashed line) and ISR/FSR (dotted line) models. 
The latter is approximately equal to the correct distribution close to the 6-quark di- 
rection and at large forward and backward rapidities. However, it lacks the t-quark 
dead-cone effect, and overestimates the radiation in the region between the b and b 
quarks. The ISR model gives a constant contribution, and is in strong disagreement 
with the correct distribution for \rj g \ ^ 2. 

Although Figs. 3 and 4 are useful are illustrating particular features of the various 
models, they correspond to idealized situations where the momenta of the t and b 
quarks are fixed. In the next section, we present more realistic distributions corre- 
sponding to the fully integrated t- and 6-quark cross sections. 

3.2 Full gluon distributions in ti events at the Tevatron 

Here we obtain the full distribution of gluons expected in ti events at the Tevatron, for 
the process qq — > ti — > bbW + W~ , with gluons generated according to Eq. (Q). We use 
MRS(H) parton distributions 0. We neglect the gg — > ti process, whose contribution 
to the cross section is an order of magnitude smaller than that for the qq initial state. 
This simplifies the interpretation of the results because the color structure of the two 
processes is different [|J. 
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Our calculations are performed entirely at the parton level, and kinematic cuts 
are kept to a minimum. We require that the b's are produced centrally and that the 
gluons fall within some typical detector pseudorapidity range. We also require that 
the gluons have sufficient transverse momentum to be detectable as soft jets but not 
so much as to invalidate the soft approximation. Finally, we require some angular 
separation between the 6's and gluons so that their respective jets are distinguishable. 
These requirements are implemented through the following kinematic cuts: 

1.5, 
3.5 , 

< 25 GeV/c , 
100 GeV , 

0.5 . (10) 

Note that the rj g cut eliminates the collinear singularities associated with initial-state 
radiation. With the exception of the cuts on AR and to some extent r) g , our results 
are not particularly sensitive to the exact values used. 

Energy and transverse momentum distributions of the resulting gluons are shown 
as solid lines in Figure 5(a) and (b). There are no surprises; the figures display the 
expected fall-off with increasing energy and px, and reflect the lower and upper cutoffs. 
Also shown in Fig. 5 are the individual contributions according to the decomposition 
in Eq. (f§), with contributions from production and decay appearing as dotted and 
dashed lines, respectively. The interference contribution is negligible for the energies 
considered here and is not shown in this plot or any that follow. We see that the 
energy distribution from the production piece extends to higher energies than that 
from decay, and that the decay pt distribution is relatively flatter than that from 
production. 

3.2.1 Angular distribution and top mass reconstruction 

The difference between the distributions of gluons produced in the production and 
decay stages is displayed more dramatically in the gluon pseudorapidity distribution, 
Figure 6. We are most interested in angular distributions in any case because we want 
to know where the soft jets are likely to appear in detectors. Figure 6 shows the net 
T) g distribution (solid line) along with the production-decay decompostion as in the 
previous figure. The decay contribution is peaked in the center and falls off quickly, 
whereas the production piece has a central dip and peaks in the forward direction. 
These results are easily understood when we consider the sources of the gluons; cf. 
Eq. (|4]). Gluons from the decay contribution come from radiation off the top and 
bottom quark lines. Radiation is largest near the quarks' direction of motion, and 
the quarks are, for the most part, produced centrally. Furthermore, we note that 



\Vb\,\Vb\ < 

\v 9 \ < 

10 GeV/c < f T 

E 9 < 

AR bg ,AR- bg > 
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the string effect in the t-b and t-b antenna pieces (the q~\pi and gjp2 terms in /'dec) 
tends to enhance central radiation compared to that in the forward regions when the 
quarks are produced centrally. Similarly, the production piece arises from radiation off 
initial-state quarks and the t and t, and we get enhancement in the forward direction. 
The net effect is a total distribution that is slightly peaked in the center and falls off 
sharply above r] g = 2 or so. 

Let us now compare these distributions with what we obtain using the simpler 
models described above. It seems plausible that the ISR/FSR model (initial-state 
and final-state radiation only, assuming color singlet initial and final states) would 
give a reasonable approximation to the correct distribution, with ISR corresponding 
to the production piece and FSR corresponding to decay. Figure 7 shows net r\ g 
distributions for the various models, with the solid line again indicating the correct 
distribution according to Eq. (|]). We see that the ISR model gives the flat rapidity 
distribution characteristic of initial-state radiation. Comparison with the production 
contribution in Fig. 6 (note the difference in vertical scales) shows that the ISR model 
significantly overestimates the amount of radiation in the central region. Adding to 
that radiation off the final 6's to obtain the ISR/FSR distribution (dashed lines) leads 
to an even greater overestimate. Hence the ISR/FSR model is not a particularly good 
approximation to the full prediction: the overall normalization is too large and the 
total shape is wrong, being too strongly peaked in the center. 

The other simple model — the BB model, in which the radiation is treated as 
in the direct bb production process qq — > bb - - does in fact reproduce the net r\ g 
distribution reasonably well; see the dot-dashed line in Fig. 7. Despite the fact that, 
for the gluon energies of interest here, the top width T t can be considered to be small, 
the T t — > oo limit seems to yield a good approximation to the full prediction, i.e., the 
results seem to be independent of the x% an d X12 factors in Eq. (|). The reason for 
this is related to the fact that radiation close to the b quarks is universal, independent 
of the process (see Ref. [Q for a discussion). But this result is deceptive for a number 
of reasons. First, depending as it does on radiation near the 6's, it is sensitive to our 
choice of AR cut; if the gluons are restricted to be farther away from the 6's, the 
agreement deteriorates. More important, although the rj distribution seems to be well 
reproduced, the distribution in azimuthal angle is not, as we shall see below. Finally, 
because radiation off the t quarks is ignored, there is no way to decompose the results 
into production and decay contributions. 

We now address the question, "What should be done with extra soft jets in attempts 
to reconstruct the top quark four-momentum?" As stated above, the decomposition 
of the gluon distributions into contributions from gluons radiated in top production 
and those radiated in decay provides guidance. Whether additional gluons in tt events 
should be combined with the b and W depends on what stage they are associated with: 
production gluons should not be combined, and decay gluons should. In principle there 
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is also a contribution from interference terms which cannot be definitely associate with 
either stage, but in practice, for gluons sufficiently energetic to be identified as soft 
jets, these interference terms are negligible. 

With this in mind, we return to the decomposition of the gluon pseudorapidity 
distribution shown in Fig. 6. We saw that gluons in the forward region are much more 
likely to have come from the production than the decay stage. Hence most forward 
gluons do not contribute to the top momentum, and based on these parton-level 
results it seems safe to conclude that soft jets with \q g \ > 1.5 can be ignored without 
introducing large systematic errors into the mass measurement. (This assumes, of 
course, that such jets are correctly identified.) Unfortunately the situation is less 
clear-cut in the central region. Despite the central suppression of production-stage 
gluons, Fig. 6 shows that gluons in the central region are almost as likely to have 
come from production as from decay; in fact, with a larger AR cut, central gluons 
would be dominated by production. Note that the situation is even worse for the 
simpler models, since the ISR rj g distribution is flat, and the BB model allows for no 
decomposition at all. 

In an attempt to identify a quantity that shows a clear excess in the decay contri- 
bution, we consider proximity to the b or b quark. It is certainly plausible to expect 
that soft jets near the fo's are likely to have come from decay-stage gluons. Figure 
8(a) shows the gluon distribution as a function of the cosine of the angle between the 
b and g, with the usual decomposition into production and decay contributions (dots 
and dashes, respectively), and total (solid). The production piece is flat, as expected 
since there is no dependence on the b direction in JF PHOD , and the decay piece increases 
as the gluon approaches the b. At first sight this looks promising - there is an excess 
in the decay piece near the b. However, this result should not be taken quite at face 
value because the amount of this excess is extremely sensitive to the cuts, especially 
that on AR between the gluon and b. The excess can be reduced or even eliminated 
by reducing the gluon range or increasing the AR cut. This would be compounded 
in practice by fragmentation of the partons into hadrons. 

We can improve the situation by taking advantage of the difference in the produc- 
tion and decay pseudorapidity distributions and making a tighter r\ g cut to eliminate 
much of the production contribution while retaining almost all of that from decay. 
Figure 8(b) shows that requiring |?7 ff | < 1.5 does indeed enhance the latter - the ex- 
cess in the region near the b quark is now more pronounced. Unfortunately, the same 
caveats about cuts and fragmentation apply here, and we cannot safely quantify this 
excess or draw any firm conclusions. 

It is clear from what we have seen above that assignments of soft jets to the 
production or decay stages cannot be made unambiguously on an event-by-event basis. 
However, we can minimize errors associated with soft jet ambiguities if we know what 
to expect from the correct gluon distribution and its decomposition and include it 
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in our simulations. Eventually we may hope to find additional discriminators that 
improve the chances of making correct assignments, leading ultimately to a useful 
prescription for dealing with soft jets. The results shown here represent a first step in 
that direction. We do not perform mass reconstruction simulations here, since we do 
not use exact kinematics in the soft gluon approximation. In future work we intend 
to dispense with the soft approximation and perform the exact calculation, which will 
allow us to investigate more specifically issues associated with mass reconstruction. 

3.2.2 Forward-backward asymmetry and color structure 

We close this section by noting that the color structure of the process process qq — > 
ti — > bbW + W~ can give rise to distinctive effects in the gluon radiation pattern 
that can be explored experimentally, at least in principle. In particular, a forward- 
backward asymmetry appears in the distribution of gluons for an appropriately chosen 
class of final-state configuations. While not directly relevant to top mass measurement, 
this asymmetry helps to illuminate further some of the physics involved, for example 
the fact that the top quarks themselves can radiate gluons. 

This asymmetry arises as a consequence of the string effects in the various antennae, 
that is, the enhancement of radiation between two quarks connected by a color string 

(see also 0). For example, in qq — > ti the qt antenna produces more radiation in 
the region between the t and q than, say, between the t and q. This gives rise to a 
forward-backward asymmetry in the gluon radiation. This asymmetry is canceled in 
the net distribution by the qt antenna, of course. However, with suitably chosen cuts 
we can enhance one contribution (from radiation off the t, say) while suppressing the 
other, thereby recovering the asymmetry. We do this by taking advantage of the fact 
that there is more radiation from a quark in the regions nearby than those far away. 

Guided also by the fact that the direction of the b quarks' momenta tend not to 
be too different from those of their parent tfs, we add to the event selection criteria 
of Eq. ( |T(]| ) the requirement that the separation in azimuthal angle between the b 
and b be greater than 135°. This tends to select events in which the parent t and 

1 have similar separation. We then preferentially select gluons associated with the t 
(as opposed to the t) by requiring that they lie within 90° in azimuth of the b quark. 
Figure 9(a) shows the resulting distribution as a solid line. A clear excess is visible 
in the forward direction, which is the direction of the initial quark's momentum. An 
equal and opposite asymmetry is obtained by considering gluons within 90° of azimuth 
of the b. 

That the source of this asymmetry is indeed the terms involving the initial quark 
and the t and i (cf. Eq. (^)) can be seen from the decomposition shown in Fig. 9(a), 
where, as in earlier figures, the production contribution is shown as the dotted line 
and the decay contribution is the dashed line. We see that the entire asymmetry 
comes from the production piece. The decay piece involves only the top and bottom 
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momenta and there is no correlation with the initial quarks. 

Finally, the asymmetry distribution reveals major discrepancies between the cor- 
rect distribution and all three of the simpler models, as shown in Figure 9(b). In the 
ISR (dotted line) and ISR/FSR (dashed line) models, in which there is no connection 
between radiation from the initial and final states, there can be no forward-backward 
asymmetry. In contrast, the BB model has a more pronounced asymmetry even than 
the correct distribution because, without radiation from the top quarks, the anten- 
nae connecting initial- and final-state quarks are enhanced. This difference between 
BB and the correct result shows that the BB model does not contain the correct az- 
imuthal dependence, despite the fact that it appears to reproduce the correct total r\ g 
distribution in Fig. 7. 

4 Conclusions 

The pattern of soft gluon radiation in tt production in hadronic collisions has a very 
rich structure. The initial-state partons, the t and b quarks can all radiate, and the 
gluon distribution in any part of phase space is a combination of what can be termed 
'production', 'decay' and 'interference' contributions. Apart from providing some in- 
teresting tests of the color structure of the events, there are important implications for 
top mass measurements from the invariant mass of the decay products. In this study 
we have attempted to illustrate the general features of the radiation pattern. We have 
quantified the idea that soft jets close to the b quarks can be associated with gluon 
emission in the top decay process, while radiation at large rapidities can be associated 
with radiation at the tt production stage. We have compared the full QCD prediction 
with those of simpler models, which are characteristic of how radiation may be im- 
plemented in Monte Carlo event generators. It is not difficult to find regions of phase 
space where the models differ significantly from the full prediction. Our conclusion is 
that the question of how to handle additional soft jets in tt events when reconstructing 
the top momentum is far from trivial. In a future study we intend to go beyond the 
soft gluon approximation to study the invariant mass distributions themselves. 
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Figure Captions 

[1] Distributions in (a) the 6-quark transverse momentum, (b) the 6-quark pseudo- 
rapidity, and (c) the azimuthal angle difference between the b and b quarks, in 
tt production, via the subprocesses qq,gg — > tt — > bbW + W~ , in pp collisions at 
y/s = 1.8 TeV. 

[2] Soft gluon radiation pattern in the rj — <p plane, for the qq — ► tt process in the tt 
center-of-mass frame. The t, i, and b, b quarks are at (r?, 0) = (0,0°), (0, 180°) 
and (r/,0) = (1,0°), (-1,180°) respectively. 

[3] The soft gluon distribution of the previous figure for fixed values of <p g . The 
curves correspond to the decomposition of Eq. (|4]): the long-dashed line is the 
'production' contribution, the short-dashed line is the 'decay' contribution, the 
dotted line is the interference contribution, and the solid line is the total. 

[4] The total soft gluon distribution of Fig. 3 (solid line) together with the dis- 
tributions corresponding to the ISR (dashed line) and ISR/FSR (dotted line) 
models. 

[5] Distributions in (a) the gluon energy and (b) the gluon transverse momentum, 
in tt production, via the subprocess qq — > tt — > bbW + W~ , in pp collisions at 
y/s = 1.8 TeV. Contributions from production (dotted lines) and decay (dashed 
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lines) are shown along with their totals (solid lines). The cuts are listed in 
Eq. ©. 

[6] Gluon pseudorapidity distributions in ti production, via the subprocess qq — > 
ti — > bbW + W~, in pp collisions at ^fs = 1.8 TeV. The net distribution is shown 
as a solid line; contributions from production (dotted lines) and decay (dashed 
lines) are also shown. The cuts are listed in Eq. fllpp. 

[7] Gluon pseudorapidity distributions in ti production, for the models discussed in 
Section 2: full distribution (solid line), ISR model (dotted line), ISR/FSR model 
(dashed line), and BB model (dash-dotted line). The cuts are listed in Eq. (|10|). 

[8] Distribution in the cosine of the angle between the gluon and the b quark, (a) 
with cuts as in Eq. ([U]) and (b) with the additional cut \r) g \ < 1.5. 

[9] Forward-backward asymmetry in gluon pseudorapidity distributions in ti pro- 
duction. The cuts are as in Eq. ([HID with the additional requirements A0 fc g > 
135° and A0j, 9 < 90°. The curves correspond to the (a) total (solid), production 
(dots) and decay (dashes) distributions, and to the distributions for the (b) full 
QCD (solid), ISR (dots), ISR/FSR (dashes) and BB (dash-dots) models. 



13 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 




QCD 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "fig2-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



This figure "fig2-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



> 

CD 

o 



tuO 

Eh 

\ 

b 
T3 



0.025 



0.020 



0.015 



0.010 



0.005 



0.000 




10 20 

P T g (GeV) 



30 



40 



0.0150 



0.0125 

CD 

o 

\ 0.0100 
^ 0.0075 

w 

\ 0.0050 



b 



0.0025 
0.0000 



20 



40 60 

E g (GeV) 



(b) 



80 



100 



Figure 5 



This figure "fig2-6.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409238vl 



0.004 



~i 1 1 r 



0.003 




a 

00 

p 

\ 
b 



0.002 



PROD 



0.001 



DEC 



o.ooo 



J I I L 



J I I L 







Figure 6 



0.006 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



0.005 



0.004 



ISR/FSR 



0.003 



0.002 



0.001 



TOT, BB 



ISR 



o.ooo 



J I I I I I I I I I I I U L 



-2 



Figure 7 



0.0030 



0.0025 



0.0020 



0.0015 



0.0010 — r 



0.0005 



0.0000 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



PROD 



TOT 



DEC 



J I I I I I I I I I I I I I I I I I L 



-0.5 



cosd 





bg 



0.5 



Figure 8(a) 



0.0030 



0.0025 



0.0020 



0.0015 



0.0010 



0.0005 



0.0000 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



DEC 



TOT 



PROD 



J I I I I I I I I I I I I I I I I I L 



-0.5 







0.5 



cos0 



bg 



Figure 8(b) 



0.0008 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



0.0006 




PROD 



0.0004 



0.0002 



DEC 



0.0000 



I - ^ fc I I I I I I I I I I I L 







Figure 9(a) 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



0.00100 



0.00075 



ISR/FSR 



BB 



TOT 



0.00050 



0.00025 



0.00000 



J I I I I I I I I I I I U L 







Figure 9(b) 



